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Abstract 

Endocytosis is an essential process by which cells internalize a piece of plasma membrane 
and material from the outside. In cells with turgor, pressure opposes membrane defor¬ 
mations, and increases the amount of force that has to be generated by the endocytic 
machinery. To determine this force, and calculate the shape of the membrane, we used 
physical theory to model an elastic surface under pressure. Accurate fits of experimental 
profiles are obtained assuming that the coated membrane is highly rigid and preferentially 
curved at the endocytic site. The forces required from the actin machinery peaks at the 
onset of deformation, indicating that once invagination has been initiated, endocytosis 
is unlikely to stall before completion. Coat proteins do not lower the initiation force 
but may affect the process by the curvature they induce. In the presence of isotropic 
curvature inducers, pulling the tip of the invagination can trigger the formation of a 
neck at the base of the invagination. Hence direct neck constriction by actin may not be 
required, while its pulling role is essential. Einally, the theory shows that anisotropic 
curvature effectors stabilize membrane invaginations, and the loss of crescent-shaped 
BAR domain proteins such as Rvsl67 could therefore trigger membrane scission. 


Author Summary 

Cells use endocytosis to intake molecules and to recycle components of their membrane. 
Even in its simplest form, endocytosis involves a large number of proteins with often 
redundant functions that are organized into a microscopic force-producing “machine”. 
Understanding how much force is needed to induce a membrane invagination is essential 
to understand how this endocytic machine may operate. We show that experimental 
membrane shapes are well described theoretically by a thin sheet elastic model including 
a difference of pressure across the membrane due to turgor. This allows us to integrate 
the different contributions that shape the membrane, and to compute the forces opposing 
membrane deformation. This calculation provides an estimate of the pulling force that 
must be generated by the endocytic machine. We also identify a membrane instability 
that could lead to vesicle budding. 


Introduction 

Endocytosis enables cells to internalize extracellular material and to recycle membrane 
components [^. During this process, the plasma membrane is deformed into an invagina¬ 
tion progressing inwards, which is severed and eventually released in the cytoplasm as a 
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Figure 1. Left : schematic representation of endocytosis. The membrane is connected mechanically 
to the actin network through sla2. Actin is polymerizing close to the basal membrane, thus generating 
an upwards (pulling) force fa. This force overcomes the rigidity hz and tension a and the turgor 
pressure 11 that all oppose invaginations. Right : four experimental membrane profiles at different 
stages, measured in an electron-microscopy study [^. Scale bar : 20nm. 


vesicle. Key endocytic components have been identified in several systems. This process 
usually involve membrane coating proteins (such as clathrin), their adaptors (such as 
epsins) and actin microfilaments together with associated factors . We focused on the 
yeast model system, in which endocytosis is well characterized experimentally. The abun¬ 
dance and localization of the principal proteinaceous components have been measured 
as a function of time both in S. pombe and S. cerevisiae [^[^. To understand how 
these components work together to deform the membrane, it is necessary to consider 
the physical constraints under which the task is performed in vivo. While in animal 
cells, invaginations are opposed mostly by membrane tension and elasticity [^, turgor 
pressure strongly opposes invaginations in plants and fungi. In those cells, the difference 
of osmolarity with the outside cause a large pressure pushing membrane against the cell 
wall. 

The magnitude of the pressure has been measured in different walled cells using 
various methods. For plant cells, studies converge to the range 0.2—1 MPa (see 
for a careful review of the subject). In the yeast S. pombe^ an effective pressure of 
0.85 ± 0.15 MPa was derived from studying the buckling of the rod-like cells in micro 
fabricated chambers [^. Based on the variation of volumes upon changes in osmolarity, 
the pressure in S. cerevisiae was recently estimated to be 0.6 ± 0.2 MPa [^, while an 
older study concluded 0.2 MPa for stationary phase cells 12 . It has been suggested that 
the pressure could be decreased locally by releasing osmolytes at the endocytic patch . 
This would however only lower the osmotic pressure, which is not opposing endocytosis. 
The hydrostatic pressure gradients equilibrate at the speed of sound (possibly 1500 m/s), 
and given the size of yeast, this is considerably faster than an endocytic event (^ 5^). 

The pressure, thus of the order of 1 MPa ^ 1 pN/nm^, pushes the plasma membrane 
outward uniformly. This effect is balanced by an equal force from the cell wall, wherever 
the membrane is in contact with the cell wall (Fig. left). The pressure strongly 
opposes endocytic membrane invagination, as membrane and cell wall must come apart. 
The large required force is produced in yeast by an actin machinery, which forms a 
crosslinked network of rigid filaments around the invagination (see and Fig. left). 
Actin polymerizes close to the basal plasma membrane [^. The newly inserted F-actin 
at the bottom of the network is though to lift the entire network away from the cell 
wall (upward on Fig. left) [^. For this force to be productive, the actin network 
should be attached to the tip of the invagination, and this is the role of the protein sla2. 
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which is required for endocytosis - while several other proteins, including clathrin, are 
dispensable [5,14. This essential role for actin is further supported by the fact that 
actin assembly precedes or coincides with membrane deformation The interplay 


between actin and pressure was nicely demonstrated by showing that impairing the 
arp2/3 complex (an actin nucleator) delayed the invagination, while decreasing the 
hydrostatic pressure (by adding sorbitol to the media) had the opposite effect 
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Coat and associated proteins play a role in endocytosis, notably by recruiting the 
endocytic machinery, and keeping the actin network physically connected to the mem¬ 
brane. Proteins that bind to the lipid bilayer can also directly induce the membrane 
to curve _ra| pr] . Two familiar examples in yeast endocytosis are clathrin and 
Rvsl67 |18|. The induced curvature is expected to be qualitatively different for these two 
proteins, because clathrin proteins form regular triskelion while Rvsl67 is shaped as 
a crescent 20 . Moreover, clathrin adaptors can bind to membrane, cargo, actin and/or 
clathrin and have essential functions in endocytosis 21 . Additionally, some adaptor 
proteins can induce membrane curvature 
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The shape of the membrane can be predicted by minimizing an effective deformation 
energy [^. This approach has been used successfully in different systems 24 26 


but the case of yeast endocytosis where the membrane detaches from the cell wall 
despite a large hydrostatic pressure has not been analyzed theoretically to our knowledge. 
Most previous work focused on the endocytosis of viruses in animal cells, which occurs 
by wrapping of the membrane around viral particles in the absence of hydrostatic 



study are consequently severely underestimated. 

In this work, we integrate the contributions of pressure, coat proteins and membrane 
properties during the endocytic invagination, providing an estimate of the force that 
actin must exert to induce the invagination. We study the effect of coat proteins such as 
clathrin that induce isotropic curvature, and contrast it with crescent-shape proteins 
such as BAR-domain that induce curvature only in one direction. Finally, we discuss 
whether the combination of actin-mediated pulling together with the removal of Rvsl67 
is sufficient to lead to vesicle internalization. 


Methods 


We predict membrane shapes by minimizing a deformation energy, using a Helfrich-type 
Hamiltonian [24p4 , in which membrane deformations are penalized by a bending rigidity 
K and tension a. We write H the difference of hydrostatic pressure between the inside 
and the outside of the cell. In addition, we assume a point force fa pulling the apex of 
the invagination, to represent the driving force generated by the actin cytoskeleton m 


We thus implicitly assume that the forces produced by actin polymerization at the base 
of the invagination, and possibly by myosin motors or other processes, are transmitted 
to the tip of the invagination over the actin network (see Fig. 1). We note hi the rigidity 
of the membrane together with its coat of proteins [l^[r^ . Moreover, we consider that 
the coat proteins curve the membrane either by scaffolding or inserting themselves in 
the membrane. We first describe proteins such as clathrin, that induce an isotropic 
curvature Co, with the same radius of curvature in both directions. 

To write the membrane deformation energy, we introduce C = ^ the local 

curvature, V the volume inside the invagination, S the surface area of the membrane 
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that is not in contact with the wall, and L the height of the invagination. The total 
energy T of an invagination reads 
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J^ = 


IL 


'liC-Cof+a 


ds + nv- faL. 


( 1 ) 


The rigidity term will tend to make invaginations as large as possible to minimize 
their curvature, while both pressure and tension will tend to make invaginations smaller, 
to minimize their volume and surface, respectively. Therefore an invagination dominated 
by pressure and rigidity will have a typical width Ru = while an invagination 

dominated by tension and rigidity will have a typical width A ^ 

From measured values of It, a (Table [^, and using the rigidity of a naked membrane 
(k ^ 40 /c^T), we find Ru ^ 3 nm and A > 9 nm. Since Ru < A in this case, we 
expect that invaginations will have a typical radius Ru and tension will not significantly 
affect the shape of the invagination. In this estimate, we have used a lower bound for 
K given by the rigidity of a pure lipid bilayer, but the actual value of k will be higher 
because it should also include the stiffness provided by coat proteins. However, the 
statement Ru < A will be all the more true for higher values of k because Ru will 
increase slower than X diS k becomes larger. We have used an upper bound for the 
membrane tension <j, but in reality it could be significantly smaller, since yeast cell have 
membrane furrows 35 , which act as membrane reservoirs and limit the surface tension. 


Moreover, we will see later how fitting the experimental membrane shapes confirms that 
the contribution of tension is negligible. 

Since the competition of pressure with membrane rigidity should control the shape 
of the invagination, the forces due to pressure and rigidity will have the same order of 
magnitude and will dominate those due to membrane tension. The total force exerted 
by pressure over the invagination should be USq, where So is the surface area of the 
wall that is not in contact with the membrane. This area is difficult to delineate in 
the electron micrographs [^, but using the dimension at the tip of the invagination Rt 
(Table to obtain an ersatz 7ri?|, we estimate the magnitude of the force to ^ 300 pN. 


39| , we can calculate 
100 nm [sl. Because 


Given this force, and a membrane viscosity r]m 10“^ kg.s-^ 
the time needed to pull a membrane tube over a distance L 
the resulting time scale Lr]/ f ^ 10“^ s is considerably shorter than the duration of an 
endocytic event (a few seconds [^), the membrane has plenty of time to reach a static 
equilibrium, at any stage of its evolution. In this quasi-static regime, we can derive the 
membrane shape equations from minimizing the energy in eq. and following earlier 
work, we also assume that they are axisymmetric for simplicity (see supplementary 
information S.I. 1.1). 

To compare our predictions to experimental membrane profiles, the measured shapes 
(Fig. a right) were projected to an axisymmetric profile (Fig. |^. For a given set of 
Co, cr, Ru^ only one value of fa allows the theoretical profile to match the invagination 


Table 1. Physical characteristics of endocytosis in yeast. 


Volume 

V - 4 X 10^ 


8 



Surface 

S' ^ 6 X 10^ nm^ 


8 



Height 

L up to 120 nm 


8 



Tip radius 

Rt ^ 12 nm 


8 



Tip curvature 

C-2/Rt-0.167 nm-i 


8 



Pressure 

H ^ 0.2—1 pN.nm~‘^ 


T 


12 

Tension 

a < 10“^ N.m~^ 


36 

37 

Rigidity 

K, > 40 kuT i.e. 160 pN.nm 




Duration 

^ b s 


s 
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Figure 2. Predicted membrane profiles (solid lines) and rectified experimental profiles (light blue 
fill) measured on electron tomograms |^. Three parameters were varied to obtain the best fit: the 
scale Ru ^ 15 — 25nm, the surface tension a and the spontaneous curvature Cq. Usually a range of 
parameter values (including a > 0) is suitable to fit each profile. The experimental membrane 
profiles (Fig. right) have been here straightened to an axisymmetric shape. Details of the method 
are found in the supplementary information (S.I. 1.5). 


shape, and therefore three parameters could be varied to fit experimental membrane 
profiles (see supplementary information S.I. 1.5 for fitting procedure). Generally, we found 
that measured profiles could be fitted precisely (Fig. except for long invaginations 
{L > 80nm, see comment in the discussion). A given experimental profiles can usually 
be fitted by a range of parameters rather than a unique set, but the values of the fitting 
parameters are remarkably consistent. The tension is negligible in ah cases as expected, 
and the values of Ru and k can be interpreted and used to calculate the magnitude of 
the apical driving force /a. We will usually express lengths as multiples of Ru and forces 
in terms of fu = dTrlli^n, since this makes our results independent of the values 11 and 
K, which are known only within an interval, without any loss of generality. 

At this stage, we emphasize that the values of k and Cq cannot be known a priori 
because these parameters represent properties of the coated membrane which is a 
complex assembly of interacting lipids and proteins, and not a simple lipid bilayer. In 
the following, we determine the rigidity n and curvature Co from the fit, assuming 
them to be constant over the membrane. In reality the coated membrane is expected 
to be inhomogeneous, but allowing these parameters to vary spatially would drastically 
increase the parameter space, such that fitting would be both underdetermined and 
computationally expensive. Moreover, the overall quality of the fits obtained by assuming 
n and Co to be homogeneous, does not warrant an increase in the complexity of the 
theory. We discuss the implications of heterogeneity in the membrane later. 


Results 

The coated membrane is rigid and preferentially curved 

The fit always provided a well determined value for Ru in the range 15—25 nm^ and 
since Ru = ^/^/2n, this allowed us to estimate the ratio between the bending rigidity 
and the pressure. Assuming 11 ~ 1 MPa one gets k ^ 2000 ksT^ while a conservative 
value of the pressure (0.2 MPa) would yield k ^ 400 kuT. The coated membrane is 
therefore much more rigid than a pure lipid bilayer {n ^ 40 ksT)^ leading to wider 
invaginations. The lower estimate corresponds to the expected rigidity of clathrin coated 
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membranes in vitro (300 ksT^ [^), while the higher value is also realistic, because other 
proteins besides clathrin are concentrated at the endocytic site 14 . Adaptor proteins 


could drastically stiffen the coat by interacting with both clathrin and the membrane . 
By increasing the effective rigidity, coat proteins enlarge the invaginations, but as a 
consequence, the overall resisting force is also increased. Importantly, our estimate of 
the minimal force needed to sustain the invagination, fa ^ 3000 pN is not based on the 
parameters that are poorly determined by the fit (such as a). 

We were also able to determine the value of the spontaneous curvature Co ^ 0.4/i?n- 
While this parameter is less well determined, it corresponds to a radius of curvature 
2/Co r\j 100 nm^ which is consistent with the known characteristics of the coat proteins. 
In particular, purified clathrin form in vitro spherical cages with an average radius of 
35 nm [^. The conditions in vivo are likely to be different however, since more proteins 
coat the membrane, and the membrane itself can also contribute to Cq. Note that we 
have assumed that curvature applied everywhere to the membrane, rather than in a 
restricted region as can be expected in reality, but this point will be discussed later. 


Large forces are required to initiate endocytosis 

By solving the membrane shape equations for a range of invagination heights L, we could 
compute the required apical force during different stages of the invagination. While 
the magnitude of the force is set by /n = dyrlli?^, its variations as a function of L are 
very instructive. The force has a non-zero value for L = 0, for the determined values of 
the spontaneous curvature Co ^ 0.4/i?n (Fig. blue curve), fa increases and reaches 
a maximum for L ^ i?n, and then decreases to a minimum value for L ^ 3i?n- For 
longer invaginations a plateau is reached. This is very different from membrane without 
spontaneous curvature (Fig. black line) or for a tension-dominated membrane 26 , 
for which the force is minimum for L = 0. In the presence of spontaneous curvature, 
the initiation force at L = 0 is high, the force peaks for smaller invaginations, and the 
pulling force is generally lowered for longer tubes (compare black and blue lines in Fig. 

i- 

The fact that an important initial force fo = fa{L = 0) is required to start the 
invagination is biologically meaningful. We estimated fo ^ 0.8/n, corresponding to 60% 
of the maximal force. Because the membrane is nearly flat if L ^ 0, we could obtain 
analytically fo = 2/n^i?n^o + 2cri?n//^ (see supplementary information S.I. 2.2). With 
(7 negligible, this expression reduces to fo = 47r/^(7o, which summarizes how coat proteins 
can affect the initial force, either by changing the spontaneous curvature, or the effective 
rigidity of the coated membrane. By itself, spontaneous curvature is not able to lift the 
membrane off the wall, but rather increases the force needed to start of the invagination. 

Because the force rises sharply with the height of the invagination, actin should be 
needed to lift the membrane already at the earliest steps of endocytosis, as recently 
observed [^[^. This idea is consistent with the observation that the delay before 
significant membrane deformation is observed depends on the competition between actin 
and pressure 15 . It has been suggested that when the membrane is still hat, actin could 


pull at one site of the membrane while simultaneously pushing on a ring-shaped zone 
surrounding this site (Fig. left). This scenario was supported by the recent observation 
that slal, an actin organizer, forms rings at endocytic sites on hat membranes, possibly 
indicating where pushing forces are applied [^. Because the membrane is still juxtaposed 
to the cell wall at this stage, the pulling and pushing forces generated by the actin 
network have to be balanced. 

Finally, a key feature of the force profile is that the most demanding part of the 
invagination occurs early during endocytosis {L < i?n)- In mechanical terms, this can 
lead to a snap-through transition 32 , in which once a critical threshold is reached. 


subsequent stages spontaneously follow because they do not require additional efforts. 
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Figure 3. Left : pulling force (normalized by /n = 27r/^/i?n) as a function of invagination length (normalized by i?n)- 
The blue line corresponds to parameters determined by fitting the experimental profiles, while the black line shows the 
behavior for an invagination without spontaneous curvature or membrane tension. Derivation of these curves is explained 
in the supplementary information. Right : membrane profiles for the force maximum (a) and for a long invagination (6). 
Lengths are given assuming Ru ^ 20nm. 


This may explain that while the duration of the endocytic early phase is highly sensitive 
to the competition between act in nucleation and hydrostatic pressure, later stages are 
largely insensitive to pressure . This snap-through force profile also elegantly explains 
why so few retraction events were observed experimentally (less than 1% of endocytic 
events fail to complete [^), despite the inherent stochasticity of a biological machine 
such as endocytosis that only contains a few hundred actin filaments [^. 

Spontaneous curvature leads to neck constriction 

Above a certain height, invagination shapes predicted by theory exhibited a neck, even 
though there is no constriction force in the model (Fig. [^. The only external force is a 
vertical lifting force applied on the apex of the invagination, and the neck appears as 
a consequence of the spontaneous curvature (Fig. last profile). Higher spontaneous 
curvatures will produce smaller necks, up to a point where the theory predicts a neck of 
zero radius, corresponding to a shape instability. 

The apparition of a localized neck as a result of global spontaneous curvature has 
been described theoretically before . We have more precisely characterized the nature 
of the shape instability (see supplementary information S.I. 3). We can understand its 
biological consequences in two limit cases : with a given value of Cq with increasing 
length L, or with a given length L and increasing Cq. 

Firstly, for small values of the spontaneous curvature, the membrane shape evolves 
smoothly upon increasing the length L of the invagination, and fa{L) is continuous. 
Above a threshold Cq however, we can distinguish two branches: one corresponding to 
short, quasi-tubular invaginations, and another one corresponding to larger spheroidal 
invaginations (Fig. |^abc, Cq = 0.45/Rn)- The two branches overlap, i.e. for a range 
of lengths, there are two possible membrane shapes for a given L, as illustrated on the 
example profiles (Fig. |^a). Above a second threshold Cq*, the two branches cease to 
overlap, such that there is a range of L where there is no equilibrium membrane shape 
(Fig. Eldef, Co = 1/Rn)- The shapes corresponding to the edges of this regions are 
illustrated in (Fig. [^d). The transition from the short invagination to the tall spheroidal 
shape corresponds to a decrease in total membrane deformation energy (Fig. Elf)- We 
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Figure 4. Illustration of the shape instability with Co = OAb/Ru (top) and Co = l/Ru (bottom). Panels a, d: profiles 
corresponding to the points indicated with the same color on the other panels of this figure. Panels b, e: pulling force 
(normalized by /n) as a function of invagination length (normalized by Ru)- With Co = OAb/Ru (top), there is an hysteresis 
with two branches of solutions: one short and tubular, the other longer and spheroidal. With Co = l/Ru (bottom), the two 
branches do not overlap and no equilibrium shape is found for a range of values of L. Dotted lines represent the unstable 
continuation of the branches. In both cases we took (7 = 0. Panels c, f: Deformation energy as a function of invagination 
length. The energy associated with membrane rigidity (i.e. bending elasticity) and the one associated with hydrostatic 
pressure are shown with discontinuous lines. The total deformation energy E = E A faL is plotted with a solid line. 


could compute the values of Cq and Cq* as a function of tension (Supplementary Fig. 4). 
The same phenomena can be observed for a given length L by increasing Cq. There is a 
threshold C^ (L) at which the predicted radius is zero, see Fig. For larger spontaneous 
curvatures Co > C^(I/), no stable membrane shape exist. This is in agreement with the 
existence of a zero-radius neck at a critical value of Co [^ . 

This shape instability is not a pearling instability (a shape instability in membrane 
tubes with tension [^), because it is inhibited by membrane tension (see supplementary 
information S.I. 3). Rather, it stems from the energetic cost of shape defects. Indeed, 
invaginations cannot be perfectly tubular, and the tip and base of the invagination 
can be viewed as defects that increase the membrane conformation energy. When the 
invagination takes an almost spherical shape, the tip defect is eliminated, and the cost 
of the base defect is minimized by making the neck infinitesimal. Indeed, we can see 
that the transition from tubular to spheroidal shapes reduces the total energy as the 
reduction in bending energy exceeds the increase of pressure-associated energy (Fig. 
f). This membrane shape instability can facilitate membrane scission, a crucial step of 
endocytosis. 


Membrane coating heterogeneity can lead to instability 


Membrane coating is not homogeneous : clathrin, sla2 locates at the tip of invagination, 
while Rvsl67 is usually located at the neck [^. Many other proteins can associate 


to membranes, including the F-Bar protein Bzzl 43 , epsins such as Entl 44 , and 
others [^|^. While we do not know the exact localization, the mechanical properties 
and the interactions of all these proteins, it is usually believed that the clathrin-coated 
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Figure 5. Theoretical membrane profiles for a constant height {L = 2.SRu) for 
increasing spontaneous curvature Co, with a = 0. 


tip is more rigid than the base of the invagination [^. This could also cause a higher 
spontaneous curvature at the tip, and additional spontaneous curvature could also stem 
from lipid asymmetry in the membrane leaflets. 

We modeled this possible heterogeneity in rigidity as a region with higher rigidity 
K (the tip), a region with a lower rigidity t^min (fhe base), and a small transition zone 
(see supplementary information S.I. 1.4). Even in the absence of spontaneous curvature, 
there is an instability if hirnin is too small compared to /^, above a threshold that depends 
on the surface area having higher rigidity. This instability resembles the curvature 
instability described above, and also stem from the destabilization of the base of the 
invagination (see supplementary information S.I. 1.4). Overall, a membrane that is more 
flexible at the base is more prone to shape instability. The results that were derived for 
a homogeneous membrane, should thus remain qualitatively valid for a heterogeneous 
membrane as well. 

BAR-domain proteins stabilize the invagination 

BAR-domain proteins are elongated crescent-like objects [^. It was shown recently that 
the number of Rvsl67 molecules, a BAR-domain coat protein, increases rapidly in the 
late stages of the endocytic process . They are thought to induce anisotropic curvature 
in the membrane. Following previous work, we assumed that they favor curvature only 
orthogonally to the symmetry axis (Z), i.e. they give a favorite radius of curvature 
Rq [^, rather and a favorite total curvature. Noting T the rigidity of a membrane 
coated with BAR-domain proteins, the contribution to the energy (eq. reads : 



This implicitly assumes that the coverage is uniform on the membrane, whereas in 
reality Rvsl67 proteins are localized at the neck of the invagination [^. This simplification 
was necessary however to derive the membrane shape equations (see supplementary 
information S.I. 1.2) and sufficient to understand in essence how anisotropic curvature 
can affect the invagination. It allowed us predict the membrane shape in the presence of 
both isotropic curvature effectors (with parameters Co, /^) and anisotropic curvature 
effectors (with parameters Rq, T)- 

We found that anisotropic curvature inhibit the membrane shape instability. The 
membrane is therefore stabilized and longer invaginations can grow continuously even 
for high values of Co (Fig. bottom, where we used Rq = Ru for simplicity). This is 
in agreement with in vivo observation that invaginations without Rvsl67 cannot grow 
longer than about 60 nm. As a corollary, removing BAR-domain proteins is a possible 
mean of triggering membrane scission. In our theory, this corresponds to T —> 0, which 
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Figure 6. Top: Theoretical membrane profiles for a membrane coated both with 
isotropic curvature effectors (such as clathrin) with Co = l/i?n and anisotropic 
membrane effectors (such as bar domain proteins) with Rq = i?n- The first three shapes 
correspond to three steps of invagination process with T = /i:. The last shape was 
obtained with T ^ 0, i.e. by removing the BAR-domain proteins. The removal leads to 
an infinitesimal neck that promote scission. Bottom: Required force as a function of 
invagination height, with and without anisotropic spontaneous curvature (T = 0 or 
using cr = 0, Co = l/Ru and Rq = Rn- la the presence of anisotropic curvature induced 
by BAR-domain proteins, there is no shape instability (the blue curve is continuous). 
Comparing with the black curve shows that the shape instability was suppressed. 


indeed destabilizes the membrane (Fig. top). This possibility is supported by the 
recent observation that membrane scission is synchronous with the disappearance of 
Rvsl67 1^. 


Discussion 

Using a general model for membrane mechanics, we could accurately fit experimental 
profiles shorter than 80nm, even though we assumed the membrane to be homogeneous, 
i.e. with constant rigidity and spontaneous curvature over the surface of the deformed 
membrane. In combination with membrane rigidity, we can expect pressure to be 
the dominant factor opposing membrane invagination during yeast endocytosis, while 
membrane tension should be negligible. This statement is derived from the dimension of 
the invagination, and from the scale of pressure determined experimentally, and is thus 
independent of the details of the model. We estimated the force required to pull the 
invagination based on the value of turgor pressure, and on the value of the rigidity that 
was determined by fitting the experimental curves. While the exact value of the force also 
depends on other parameters and in particular on the height of the invagination (see Fig. 
1^, its scale is primarily determined by pressure and the width of the invagination. For 
the measured range of pressure 11 ^ 0.2—1 MPa, the force scale is /n ^ 1000—5000 pN. 
This is significantly larger than previously estimated (1 — 1000 pN 13,48,^). The 
corresponding range of value for the rigidity is ^ 400—2000 ksT. This is much stiffer 
than a pure lipid bilayer, because the membrane is heavily coated at the endocytic 
site. It was suggested that phase boundaries could play a role by generating a line 
tension, favorable to membrane budding. However the typical line tension (of the order 
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of ^ 0.4 pN [^) is much smaller than /n, and such phenomenon can therefore be 
discarded. 

The fact that the scale of the invagination is determined by the ratio of rigidity and 
pressure {Ru = ^k,/2IV) suggests a possible tradeoff. In the absence of any reinforcement, 
the naked membrane could only make small invaginations with a width of 3 nm. The 
coat enlarges the invagination by a factor ^ 10 by increasing the local rigidity, but at 
the same time also increases the required pulling force. The mechanical properties of 
the coat thus likely represent an optimal tradeoff between increasing the radius of the 
invagination (making individual endocytic events more productive) and limiting the 
driving force required from the actin machinery. The ratio k/YL is likely to have been 
adjusted in the course of evolution, and is not expected to be conserved across species. 


In general, coat protein can induce a negative tension that can favor tabulation 51 
because of their adhesion energy with the membrane. This adhesion energy has been 
estimated to cu ^ 10~^N/m in the case of clathrin 52 , while other coat protein may 
have adhesion energies of uj r\j 10 ^N/m Izl Using the largest of the two values, we 
find a pulling force ttujR ^ 30pN^ which is insufficient to drive the invagination against 
the turgor pressure. Moreover, we found that spontaneous curvature actually increases 
the initial pulling force in these conditions, and consequently coat proteins do not 
help to initially lift the membrane. This counter-intuitive result is in agreement with 
observations that clathrin is not necessary to initiate curvature during endocytosis in 
yeast [^. This result is specific to the situation where a membrane subjected to a 
large pressure has to be pulled away from its supporting wall. In this configuration, 
the invagination must have regions with positive (at the tip) and negative (at the rim) 
curvature, and the energetic cost of the latter are significant in the presence of positive 
spontaneous curvature. 

The good quality of the fits (Fig. indicates that the experimental invagination 
profiles are consistent with a pulling point force, directed away from the cell wall. In the 
cell, the forces generated in the actin network are probably transferred to the tip of the 
invagination over an finite area corresponding to the sla2 proteins. In the model, the 
point force is effectively propagated over the tip of the invagination, due to the inferred 
rigidity of the coated membrane. As long as the coat structure remains sufficiently rigid, 
these two situations should be equivalent. 

We find that under physiological conditions of pressure, direct constriction by actin 
in a mechanism similar to cytokinesis is not necessary to explain the shape of the 
invagination, or specifically the formation of a neck. This result is similar to what was 
reported in the absence of turgor 29 , and stems from having an isotropic spontaneous 


curvature that is everywhere the same on the membrane. 

With isotropic spontaneous curvature, the membrane can have a shape instability 
(Figs. 4 and 5). A more complicated model that included some level of inhomogeneity 
(S.I. 1.4), showed that this shape instability is facilitated if the coated membrane is 
less rigid at the base than at the tip of the invagination (S.I. 3.2). We thus expect the 
shape instability to be promoted by the inhomogeneity of the system, in particular the 
fact that the membrane tip is covered with clathrin while the base most likely is not. 
Qualitatively, pressure is pushing on the neck, thus constricting it, while membrane 
tension is pulling on it. Thus the shape instability is promoted by pressure and inhibited 
by membrane tension (S.I. 3.1). 

Crescent-shaped curvature effectors such as Rvsl67 were modelled as inducing 
anisotropic curvature. The analysis showed that their presence stabilizes the invagination, 
suggesting that Rvsl67 removal can cause membrane destabilization and scission. It is 
known that membrane fission is promoted by amphipathic helix insertion and inhibited 
by crescent BAR domains 53 , and a theoretical argument was made comparing idealized 
tubes to spherical vesicles. Here we show that crescent-shaped proteins also stabilize 


11/23 










intermediate shapes of endocytosis, in what appears to be a very generic phenomenon. 

It was reported that Rvsl67-deprived cells have shorter invaginations [^. A possible 
explanation is that BAR domain proteins stabilize invaginations, as in our theory. An 
alternative possibility is that Rvsl67 would be needed to lower the pulling force required 
at later stages of endocytosis. However, our calculations indicate that the highest forces 
are required during the initial stages of endocytosis (Fig. |^. Indeed experimentally, the 


growth speed of the later stages seems independent from the pressure 15 , which is the 


main force opposing endocytosis. This indicates that stabilization of the shape by BAR 
domains, rather than reduction of force, is the most likely explanation of the observed 
shortening of invaginations in Rvsl67 mutants. 

The homogeneous theory could not fit large membrane deformations (> 80nm). 
Indeed, Rvsl67 appears at the neck for large invaginations [^, and probably induce 
inhomogeneous properties near the base. This effect was not included in the inhomoge¬ 
neous model described in (S.I. 1.4), because introducing several different regions makes 
it technically intractable. 

Endocytosis is ubiquitous in nature and in other model systems such as animal 
cells, it appears to involve a similar set of molecules. The physical conditions can 
however vary greatly in different cells, and this should be reflected in some of the 
requirements that have constrained the evolution of the endocytic machinery. In the 
absence of significant hydrostatic pressure, membrane coating is sufficient to generate 
invaginations [^, while dissipative processes involving actin and/or dynamin machinery 


appear necessary to membrane severing 54 55 . For the yeast system that we examined. 


pulling the membrane away from the cell wall seems sufficient to induce a complete 
budding event. The next crucial step is to decipher the mechanism by which the actin 
cytoskeleton is able to exert the required force. Precise quantitative information is 
available to do so and theoretical work is under way [^. Ultimately, this will 
offer a unified model incorporating both actin cytoskeletal dynamics and membrane 
mechanics, based on the experimental observations. 
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1 S.I. : Membrane shape 

1.1 Differential equations 

The membrane shape equations correspond to the minimum of the membrane energy 
(main text, equation 1). To derive these equations, we renormalize energies by n 
and distances by without loss of generality. The typical scale of forces will be 


/n = 27r/^/i?n- Following earlier work 34 , we assume rotational symmetry around the 
vertical axis, which is defined as the direction normal to the cell wall, orientated towards 
the cell inside. Let Z be the distance from the cell wall along this axis, R the distance 
from the symmetry axis, s the arclength along the membrane, and ip the angle of the 
membrane with the horizontal axis (see Fig. left), all being dimensionless quantities. 
Additionally, we define a = Ru^af f = fa/ Co = RuCq and H = 1/2. To impose 
the geometric relations dR/ds = R = cos {pj) and Z = — sin (t/), we use two Lagrange 
multipliers u and rj. The renormalized energy R' finally reads : 


R' = 27r f C{R,2p,ip)ds, 

Js 


(3) 


with C the renormalized energy density: 


R ( sin lb - \ ^ YIR^ 

jC= — lip^ - R -+ + ’ ( 4 ) 

where the dot indicate the derivative with respect to the arc length s. From these, we 


derive the Euler-Lagrange equations 34 


cos 7 / / sin 7 / ; \ n ^ ^ v , ^ 

ip = —^ - 'ip + V- f j ^ -Rcosip^ -^smip, 

^ = pi’ - Cof - +a + nRsm'tp, 

ri = 0, 


(5) 

(6) 
(7) 


which are valid between R = 0 (the center of the invagination) and R = Ri^ after which 
the membrane is in contact with the cell wall. Since the energy involves ip^ ip has to be 
continuous and hence ip = 0 sit R = Ri^ whereas ip can be discontinuous. Canceling the 


boundary term 


Hamiltonian 


34 


uRi] imposes z/(0) = 0 (as Ri is a free parameter) and canceling the 


imposes ip{0) = — 


Cq + 2(7 . 


For simplicity, we assumed Co to be 


constant over the whole membrane surface, an approximation discussed later. 

While assuming axisymmetry may seem a strong approximation, minimizing Helfrich’s 
Hamiltonian will produce axisymmetric shapes, and a more general theory would yield 
the same result, unless the force fa is not orthogonal to the cell wall, or if the membrane 
properties are severely anisotropic [^. Experimentally, we see membrane profiles that 
are not perfectly axisymmetric [^, which is expected considering the small size of the 
actin machinery, where possibly as few as 10 filaments are polymerising at any time |^. 


1.2 Membrane with anisotropic coat 

To describe a membrane coated with anisotropic coat proteins such as Rvsl67, we added 


a term favoring a radius Rq to the energy density 47 




(8) 
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Cell wall 


Figure 7. Diagram of the endocytic process and its parametrization. We assumed 
rotational symmetry around the Z axis, normal to the cell wall. The shape is 
parametrized by the curvilinear abscissa 5, the distance R from the symmetry axis, and 
the angle ip. L is the maximal height of the invagination, Ri is the distance beyond 
which the membrane is in contact with the cell wall, and 11 is the difference of 
hydrostatic pressure across the plasma membrane. 


in which is the new renormalized energy density, T = r/K and Rq = Rq/Ru • 

The differential equation for ip is unchanged (equationbut a term ^T{l/R — 1 /Rq)^ 
is added to z> (equation . Additionally, the boundary value of ip is now : 

^(0) = -Jc7oV2<7 + f ^^~y (9) 

V R^Ro 

1.3 Numerical solutions 

To find the membrane shape, we integrated equations [5][7| from R = Ri down to ^ 0, 
i.e. taking a negative ds. Because of the divergence at i? = 0, in practice we integrated 
from Ri to a non-zero cutoff e = 0.004, that was chosen to be smaller than a protein. 
We used a shooting method to find the values of Ri and / such that the membrane 
reaches Z = L and ip = 0 ai R = e. Because of the strong non-linearity in the equations, 
most initial values of Ri and / will lead to a diverging shape. However, by scanning 
randomly for initial values of Ri and / for a given L, we could find parameters that gave 
a physical solution to the equations. 

To compute membrane shape for a range of L for a given value of Co and d, we started 
from a small value of L and determined the best {Ri, /} by a random search coupled to a 
shooting method. Then, we looked for the solution for L-i-SL using the shooting method 
starting from the same Ri and /. When SL was small (typically 10“^—the 
shooting method usually converged, unless a bifurcation was reached. To find solutions 
for other parameters Cq and d, we either started another random parameter search 
coupled to the shooting method, or started from a known solution, and slowly changed 
the Co and/or d (again with small steps). 

1.4 Heterogeneous membrane 

To model the existence of membrane heterogeneities, namely the clathrin-rich tip of the 
invagination, we introduced a non-constant rigidity k i-d Ka{s). Assuming a{s) to be 
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Figure 8. Top : three membrane shape profiles. Thin line: Co = 0,a = 0, dashed line: Co,d = 1.0, and thick solid line: 
(7o,d = 0. Bottom : force versus height of the invagination, for the three aforementioned parameter sets. The gray and noisy 
lines represent the numerical calculation of dE/dL^ showing that our results are thermodynamically coherent. Dotted lines 
indicate the analytical predictions of the plateau force. 


derivable, this changes the differential equations for t/; : 


^ = 


cos^ip f sin?/; 


R 


R 


- + 


V- f 


a 


n 


a 


+ —Rcos^p H- — sin?/;-+ 


2a 


aR 


a 


sin?/; 


-Co (10) 


^ - Co)^ - + cr + nRsmi) (11) 

f] = 0. ( 12 ) 

For (a(s) we used a function varying smoothly between 1 and over a width As : 

(a(s) = 1 if s < Si (13) 

if s > Si (14) 


a{s) = 1 -h (a^ 


. , (si-sy 

1) ( 1 — e 2 As 2 


Here, si is the position of the step along the arclength (with s = 0 dit R = Ri and 
s = Stot Sit R = e). We determined si numerically to obtain that the surface area defined 
by s < Si is the tip surface Sup. 


1.5 Fitting of experimental data 

The 2D experimental profiles derived from electron microscopy imaging are available 
online. They needed to be projected to an axisymmetric profile, since the theory could 
only consider axisymmetric shapes. For a given profile, we find the two points A,B 
where the membrane intersects an horizontal axis, at each height from the base. We 
then extract the two angles of the membrane with respect to the vertical axis at points 
A and B, and we define the local orientation t as the average of these angles. The radius 
of an axisymmetric profile at this height is then [AB] cos (t)/2. 

In the fitting procedure, it is important to note that Ru act as a scaling parameter, 
while a and Cq control the shape of the invagination. The renormalized force / can 
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be seen as a Lagrange multiplier controlling the height of the invagination. Therefore, 
we compared experimental shapes to theoretical profiles with a large range for d and 
Co, and for each profile, we tried values of Ru from 5nm to 35nm. For a given set 
{Co,d,i?n }5 only one value of fa produces the correct invagination height. 

The error between a theoretical profile and an experimental shape was defined as : 


Error 


(15) 


Where the Ri are the radii of profile sections, for height between 0 and L, where L is 
the length of the invagination. 

A range of parameters {Co, d, i^n}, rather than a discrete set, can be used to fit each 
profile, while keeping the error below 5%. We found that most profiles could be fitted 
with Ru ^ 15 — 25nm, d ^ 0, and Co ^ 0-0.5. The accuracy of the fit is undermined 
by the facts that we did not fit the original profiles but an axisymmetric projection. 
Moreover the solution space is degenerate, and bias are possible in the experimental 
methods. 


2 S.I. : Forces 


2.1 Force as the derivative of the energy 


In our physical theory, the force fa pulling the invagination should satisfy the thermody¬ 
namic definition of a force, i.e. fa = —dE/dZ, in which E = EfaL is the total internal 
energy. Hence / can also be seen as a Lagrange multiplier controlling the height of the 
invagination. We have several ways to test the validity of our force computation. First, 
we can compute analytically the force required to pull long invaginations {L Ru): 
that are mostly tubular [^. In this configuration, the total energy of the system is 
dominated by the energy of the tubular portion, and the radius of the tube i?* can be 
found by minimizing the energy per unit length of tube: 


Efube — ‘^TT R 




HTTi?^ 


^R^tube 


= 0 . 


(16) 

(17) 


R* 


The numerically computed plateau value for the force matches this analytic value 
/* = Etube{R*) perfectly (figurebottom). 

Second, to confirm that the force computation is exact for any L, we verified that the 
input value for fa, found by the shooting method, corresponded to the thermodynamical 
definition of the force OlE. The numerical estimation of d^E/fu is noisy (due to small 
error of on the height L and on the angle at the tip), but indeed / always matches 
OlE/ fu (figure]^ right, underlying grey lines). 


2.2 Initiation force 

At the early stages of endocytosis, before the appearance of a neck, we can describe the 
membrane in the Monge form using Z{R) and Z' = dZjdR [^. For early invaginations, 
we can further take the small deformation limit Z' ^ 0, leading to a simpler formulation 
of the energy: 


E 


M 


27TR 


pRi - 

L 


(AZ - Cof2aZ'‘^ + 2nZ 


RdR . 


(18) 
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Figure 9. Initiation force /o as a function of physical parameters Co and d. Black line 
: numerical solution in the Monge approximation. Crosses : numerical solution of the 

general shape equations. Grey line : fit by /o = 2/n A/ Cq^ + 2d. 


Taking the limit a 
length L : 


0, equations 5|6 become, for an invagination of renormalized 


Z{R) = -^R^{R^+ ^2a)+L + hR^\og{R). (19) 

izo 

The parameters {a, 6, should satisfy the boundary conditions : 

Z{Ri) = Z'{Ri)=0, (20) 

Z"{Ri) = -\/Co^ + 2a. (21) 

For any given Z, we can then compute numerically Z{R)^ from which we deduce Em{L). 
We can then compute the force OEm / dL numerically. We could compare the approximate 
result from the Monge representation (in the absence of tension), the numerical result 
from the general membrane shape, and we find that they are well fitted by the formula : 

/o = 2/n7(7o"+2d (22) 


3 S.I : Shape instability 


3.1 Curvature-induced shape instability 

It is interesting to compare the shape instability that we found with the pearling insta¬ 


bilities (see 59 and references therein), whereby (local) spontaneous curvature competes 
with membrane tension , usually under constraints of volume conservation . 

Comparing the energy of a tube (equation with Z = '7r/2 and Z = 0) to that of a 
string of spheres does not yield any instability in our case. Instead, upon increasing L 
at the transition, the bending energy decreases more than the pressure energy increases, 
and this transition occurs at heights where the membrane is far from tubular. This 
transition is therefore caused by the cost of bending the membrane at the tip and tail of 
the invagination, and these defects cost more energy in the quasi-tubular shape than in 
the spheroid shape, where the neck is a defect that can become infinitesimally small. 
The phase diagram (Fig. right), shows that membrane tension actually inhibits this 
transition, which is due to membrane tubes of minimum energy having less membrane 
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Figure 10. Left : Force as a function of spontaneous curvature for L = 2.SRu and a = 0; the close-up shows that the 
branch ends at a finite Cq. Right: Phase diagram of the budding instability as a function of spontaneous curvature Cq and 
tension a. Grey: no instability, green: instability with hysteresis, red: instability with discontinuity (lack of solutions for a 
range of L). Solid lines are visual guidelines showing that the critical value of Cq depends linearly on a. 


area than spheroids. The phase diagram also exhibit a region in which there is no 
stable membrane shape for a range of invagination heights, in contrast with the pearling 
instability. Despite the complexity of the equations, the phase boundaries exhibit a 
linear behavior for the critical value of Cq as a function of a (Fig. 


10, right). 


3.2 Rigidity-induced shape instability 

Strikingly, there is a shape instability even in the absence of spontaneous curvature if 
the membrane is heterogeneous, with a basal membrane of lesser rigidity 
This instability appears if nmin !is smaller than a certain threshold, which depends on 
the surface area with higher rigidity. For example, with tip surface of we found a 

threshold of 0.25/1:. 

We then used K.rnin = O-l/i: (representative of the ratio of clathrin rigidity to membrane 
rigidity). Though the membrane shapes are greatly influenced by membrane heterogeneity 
(Fig. right), their behavior as well as that of the force-distance curve (Fig. left) 
is very similar to the curvature-induced instability, and it appears that both instabilities 
are of the same nature. 
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